The Space Density of Primordial Gas Clouds near Galaxies and 
Groups and their Relation to Galactic HVCs 
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ABSTRACT 

The Arecibo H i Strip Survey probed the halos of ~300 cataloged galaxies and the envi- 
ronments of ^lA groups with sensitivity to neutral hydrogen masses > 10^ Mq. The survey 
detected no objects with properties resembling the High Velocity Clouds (HVCs) associated with 
the Milky Way or Local Group. If the HVCs were typically Mhi = 10^'^ Mq objects distributed 
throughout groups and galaxy halos at distances of '^1 Mpc, the survey should have made '^70 
HVC detections in groups and ^^250 detections around galaxies. The null detection implies that 
HVCs are deployed at typical distances of <200 kpc from the galaxies or group barycenters. If 
the clouds are in virial equilibrium, their average dark matter fraction must be 98% or higher. 



Subject headings: 
mass function — 

1. Introduction 



Local Group — intergalactic medium - 
radio lines; ISM 



After an extensive review of the High Veloc- 
ity Cloud (HVC) literature, Wakker &; van Woer- 
den (1997) concluded that no single origin can ac- 
count for the properties of the HVC population of 
neutral hydrogen clouds that surround the Milky 
Way Galaxy. Instead, several mechanisms, includ- 
ing infalling extragalactic clouds, cloud circulation 
within the Galactic halo driven by a galactic foun- 
tain, and a warped outer arm extension to the 
Galaxy must be invoked. The Magellanic Stream 
and associated HVC complexes require an expla- 
nation due to tidal interactions within the Local 
Group (LG) (cf. Putman et al. 1998). 

A defining property for the HVCs has been 
the lack of associated stellar emission. This also 
means that spectroscopic parallax methods can- 
not be used to measure distances to the clouds, 
and the lack of known distances, in turn, hinders 
the calculation of the clouds' physical properties. 
Recently, absorption lines (cf. van Woerden et 
al. 1999) and Balmer recombination line emission 
driven by reprocessing of the ionizing radiation 
originating in the Galactic star forming regions 
(Bland-Hawthorn et al. 1998; Bland-Hawthorn & 
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Maloney 1999) have been used to specify distances 
to a few of the clouds, indicating a range of dis- 
tances from within the Galactic halo to greater 
than 50 kpc. 

Recent interest has returned to the idea that 
HVCs could be primordial objects raining on the 
Galaxy, as either remnants from the formation of 
the LG or as representatives from an intergalactic 
population of dark matter dominated mini-halos 
in which hydrogen has collected and remained sta- 
ble on cosmological time scales (Blitz et al. 1999, 
BSTHB; Braun & Burton 1999, BB). The require- 
ment of gravitational stability without star for- 
mation places lower limits on the distances of the 
clouds from the Sun - a sort of independent dis- 
tance indicator that can be applied to each cloud 
individually, depending on its H i flux, angular 
extent and velocity profile width. Typical derived 
distances are of order 1 Mpc, implying that this 
segment of the HVC population (1) inhabits the 
LG rather than the Galactic halo, (2) has typical 
H I mass per cloud greater than ~10^ Mq, and (3) 
increases the LG H i mass budget by contributing 
-4x10^° Mq in the case of the BSTHB sample. 
For the BB sample, which is restricted to 65 con- 
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firmed "compact HVCs" , the integral H i content 
adds ~109 Mq to the LG. 

Further impetus to search for a primordial pop- 
ulation of low mass objects comes from simula- 
tions of galaxy and group formation (cf. Klypin 
et al. 1999), which predict of order 10 times 
more ~10^ M© mini-halos in the LG than can be 
counted in the dwarf galaxy population. The as- 
sociation of HVCs with this missing population, 
as well as arguments based on the kinematics of 
the cloud population (BSTHB and BB), makes an 
appealing picture for the extragalactic/LG expla- 
nation. 

Similar concerns motivated the extragalactic 
21cm line survey by Weinberg et al. (1991), whose 
study of several representative environments (clus- 
ters and voids) found only gas-rich galaxies con- 
taining stars. Several extragalactic H i surveys of 
substantially larger volumes to more sensitive H i 
mass limits have also found no objects with HVC 
properties (i.e. H i detections without associated 
starlight) (Zwaan et al. 1997; Spitzak & Schneider 
1999; Kilborn ct al. 1999). 

Clearly, if the HVC phenomenon is a common 
feature of galaxy formation and evolution, then 
extragalactic surveys of the halos and group envi- 
ronments of nearby galaxies should show evidence 
for this population. We take two approaches to the 
problem of placing the local HVC population in an 
extragalactic context. The first (in sections 3 and 
4) is to compute the H i mass function for the local 
group, both for optically selected group members 
and for group members plus HVC populations as 
modeled by BSTHB and BB. The second, separate 
approach (section 4) is to calculate the probability 
that the narrow strip that the Arecibo^ H i Strip 
Survey (AHiSS) (Sorar 1994; Zwaan et al. 1997) 
makes through the halos of ^200 galaxy halos and 
~14 groups would detect members of HVC popu- 
lations in those systems. 

2. The Local Group H i Mass Function 

The H I mass function (HiMF) is used in extra- 
galactic astronomy to quantify the space density of 
gas-rich galaxies and possible intergalactic clouds 

The Arecibo Observatory is part of the National Astronomy 
and Ionosphere Center, which is operated by Cornell Uni- 
versity under a cooperative agreement with the National 
Science Foundation 



as a function of H i mass. For the field galaxy pop- 
ulation, the HiMF has been determined accurately 
for Mm > 10^-^ h^^ Mq and can be fit satisfac- 
torily with a faint end slope a « —1.2 (Briggs & 
Rao 1993; Zwaan et al. 1997; Kilborn et al. 1999). 
At lower Mhi < 10^'^ h^^ M©, where there is con- 
siderable uncertainty due to the small number of 
detections, Schneider et al. (1998) have found ev- 
idence for a steep upturn in the tail of the HiMF. 
Although this steep tail has a tantalizing similar- 
ity to the signature of massive HVCs, it appears 
that at least one of the two H i signals responsible 
for the rise comes from a normal galaxy, and the 
other is too close to a bright star to exclude faint 
optical emission (Spitzak & Schneider 1999). 

We construct the HiMF for the LG from H I 
measurements of all known LG members as com- 
piled recently by Mateo (1998). Within the LG 
volume, a total of 22 galaxies arc known in which 
H I has been detected. The statistics are therefore 
poor, and consequently the HiMF is noisy. The 
LG HiMF is in one sense the best measured HiMF, 
with data over six orders of magnitude, compared 
to three or four orders of magnitude for determina- 
tion for the field galaxy population. On the other 
hand, the LG HiMF may suffer from severe selec- 
tion effects due to obscuration by dust in the disk 
of the Milky Way galaxy. In order to estimate 
how many galaxies might have escaped detection 
so far, Mateo (1998) plotted the cumulative num- 
ber of galaxies as a function of 1 — sin for the 
total sample of LG galaxies. If LG galaxies were 
distributed equally over the sky, the resulting his- 
togram should be a straight lino. Wo applied the 
same method to the sample of galaxies with H I 
detections and found that 4 to 7 galaxies with H i 
masses > 10^ Mq arc likely to be missing at low 
Galactic latitude. The missing galaxies would be 
predominantly the ones with low optical surface 
brightness, but since there is no clear correlation 
between surface brightness and H i mass, it is not 
possible to make a more refined correction to the 
HiMF than just adding one galaxy to each half 
decade mass bin below Mhi = 10^ M© . For larger 
galaxies with Mhi > 10^ M©, we assume that the 
census of the LG galaxy population is complete 
(cf. Henning et al. 1998). 

The result is presented in the top panel of 
Fig. 1, where the points represent the LG HiMF 
and the line is the field HiMF derived by Zwaan et 
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Fig. 1. — Top panel: H I mass function of the Local Group 
(LG). The points show the space density of LG members 
containing H I, after correcting for incompleteness. The 
solid line shows the field HiMF from Zwaan ct al. (1997), 
scaled vertically so as to fit the points. The dotted line 
is a HiMF with a steep upturn at the low mass end, re- 
cently proposed by Schneider et al. (1998) Second panel: 
H I mass functions for extragalactic HVCs. The dark his- 
togram shows the space density of BSTHB HVCs if they 
are put at the critical radii for gravitational stability. The 
unshaded histogram shows the effect of a baryon to to- 
tal mass ratio / varying from cloud to cloud. The lighter 
shaded histogram corresponds to BB HVCs, all at the same 
distance of 1 Mpc. Lines are as in top panel. Third panel: 
H I mass function for BSTHB HVCs for different values / 
(from right to left: /=0.2, 0.1, 0.05, 0.025, 0.0125). The 
HVC HiMF is consistent with the field HiMF if / < 0.02 
and the median distance < 200 kpc. Bottom panel: Distri- 
bution of the expected H l detections in the AHiSS for the 
BSTHB population. The light grey histogram corresponds 
to clouds surrounding galaxies, the dark grey histogram 
represents clouds in galaxy groups. See section 4 for expla- 
nation. 



al. (1997) scaled vertically so as to fit the points 
in the region around Mhi = 10^ Mq where the 
curve has been measured accurately. This scaling 
is justified given the fact that the HiMF for opti- 
cally selected and H i selected galaxies is identi- 
cal (cf. Briggs & Rao 1993; Zwaan et al. 1997). 
The scaling accotmts for the ovcrdcnsity of the 
LG, which amounts to a factor of 25, assuming 
that the LG volume is 15 Mpc^. Also shown is a 
HiMF with a steep upturn proposed by Schneider 
et al. (1998). The large divergence between the 
extrapolated curves from Zwaan et al. (1997) and 
Schneider ct al. (1998) illustrates the uncertainty 
in HiMF below Mm = lO^-^ Mq. The LG HiMF 
for optically selected galaxies is remarkably flat, 
with the faint-end slope of a Schechter function fit 
of a « -1.0. 

Studies of the H l content of galaxies in dif- 
ferent environments (including voids [Szomoru et 
al. 1996], clusters [McMahon 1993] and groups 
[Kraan-Korteweg et al. 1999]) have shown that the 
shape of the HiMF for Mhi > 10® Mq is indepen- 
dent of cosmic density. The fact that we find here 
that the HiMF of optically selected LG members 
is flat down to H I masses of a few x 10^ Mq does 
not, however, insure that the field HiMF is also fiat 
to these low masses. Although the crossing time of 
the LG is approximately equal to a Hubble time, 
there are clear indications of interactions (Mateo 
1998, and references therein). The H i distribu- 
tions in the lowest luminosity LG dwarfs are often 
highly asymmetric, perhaps indicative of tidal dis- 
tortions. It is quite possible that low mass systems 
arc destroyed or merged, which could cause the LG 
HiMF to be fiatter than that of the field. 

3. H I Mass Functions for Extragalactic 
HVCs 

We now derive HiMFs for the population of ex- 
tragalactic HVCs as proposed by BSTHB and BB. 
For the compact, isolated HVCs identified by BB, 
we estimate H I masses by using their measure- 
ments of the integrated fluxes and the assump- 
tion that all clouds are at the same 1 Mpc dis- 
tance, suggested by BB. Typical H i masses are 
then a few times 10^ Mq and the sizes are ap- 
proximately 15 kpc. The resulting HiMF is indi- 
cated by the light shaded histogram in the second 
panel of Fig. 1. For the BSTHB sample we use the 
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compilation of HVCs by Wakker & van Woerden 
(1991, WW), which forms the main source for the 
BSTHB analysis. For each cloud the distance rg at 
which the cloud is gravitationally stable is calcu- 
lated. Adopting BHTHB's value of / = 0.1 for the 
ratio of baryonic mass to total mass leads to typi- 
cal distances of self-gravitating HVCs of ~ 1 Mpc. 

The dark shaded histogram represents the 
HiMF for the BSTHB clouds, assuming that 
all clouds are at their distance Vg. The result- 
ing HiMF shows a clear peak at approximately 
3 X 10^ Mq, equal to the typical H i mass esti- 
mated by BSTHB. The largest iuic;crtainty in the 
determination of Mhi is /, which might vary from 
cloud to cloud. We tested the effect of a varying 
/ on the HiMF by applying to each individual 
cloud a random value of / in the range 0.03 to 
0.3. The resulting HiMF (shown by the unshaded 
histogram) is not significantly different from the 
HiMF based on a fixed / in the region of interest 
> IO^Mq. 

The space density of the BB HVC population is 

comparable to the HiMF for optically selected LG 
galaxies, but it is obvious that the BSTHB HVCs 
outnumber normal galaxies by a factor 5 to 10 in 
the range 10^-^ < Mhi < lO'' Mq. This implies 
that if the BSTHB cloud population is typical for 
galaxy groups, H I surveys in groups should have 
encountered 5 to 10 dark H i clouds for every de- 
tected galaxy. This is clearly at variance with the 
observations. 

At what distances from the Milky Way must 
the HVCs be located so that their HiMF is not 
in conflict with the observed field HiMF? Since 
the virial distance rg is directly proportional to /, 
we can test this by varying / from 0.2 to 0.025 
and calculate the resulting HiMF. The results are 
shown in the third panel of Fig. 1. The space 
density of HVCs in the LG can only be brought 
into agreement with the observed field HiMF if 
the median value of / is lowered to ^ 0.02, a 
value much lower than what is normally observed 
in galaxies. The median distance of such clouds 
must be smaller than ^ 200 kpc. 

4. Expected number of extragalactic HVC 
detections 

The Arecibo H i Strip Survey (AHiSS), which 
is discussed in detail in Sorar (1994) and Zwaan et 



al. (1997), puts limits on the space density of pri- 
mordial gas clouds in external galaxy groups and 
around galaxies. The survey was taken in drift- 
scan mode and consists of two strips at constant 
declinations, together covering 20*^ of RA over a 
redshift range from to 7500 kms~^. The limit- 
ing column density was 10^^ cm~^, which is lower 
than that of most of the HVCs in the WW com- 
pilation and those presented in BB. The sky cov- 
erage is small (15 deg^ excluding the side lobes) 
but the survey strip passes through the halo re- 
gions of many groups and galaxies as shown in 
Fig. 2. The unique character of this survey makes 
it therefore more suitable to assess the HVC prob- 
lem than other surveys of equal size. Prom the 
Lyon-Mcudon Extragalactic Database (LEDA^) 
we selected all known galaxies with projected dis- 
tances < l/igs^ Mpc from the two AHiSS strips. 
The circles in Fig. 2 indicate shells with 1 Mpc 
radii around the galaxies. Since the discussion in 
BSTHB is primarily focussed on galaxy groups, we 
also selected all cataloged groups within 1 Mpc of 
the strips. Galaxy groups were drawn from Willick 
et al. (1997), who used the the Mark III catalog, 
and Garcia (1993) who selected groups from the 
LEDA galaxy sample. 

We fill the volumes around the selected groups 
and galaxies with a synthetic population of HVCs 
similar to the one proposed in BSTHB. To con- 
struct such an ensemble we make use of the com- 
pilation of HVCs by WW as discussed in sec- 
tion 3. Although BSTHB put particular emphasis 
on galaxy groups, we choose to consider clouds 
around individual galaxies as well. Hierarchical 
formation scenarios do not distinguish between 
galaxies and groups in the relative number of satel- 
lites (Klypin et al. 1999). Further motivation 
comes from the fact that in the LG, subcluster- 
ing of dwarfs is observed around the Milky Way 
and M31 (Mateo 1998). 

The 3' beam of the Arecibo telescope subtends 
c^beam — 0.871? kpc at a distance D Mpc. The 
typical sizes of the HVCs discussed in BSTHB are 
28 kpc and the lowest column density clouds could 
therefore be detected out to distances of 32 Mpc, 
beyond which the average HVC would no longer 

^We have made use of the Lyon-Meudon Extragalactic 
Database (LEDA) supplied by the LEDA team at the 
CRAL-Observatoire de Lyon (Prance). 
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Fig. 2. Illustration of the two surveys strips and 
the distribution of spheres with radii of 1 Mpc around 
known galaxies (thin circles) and galaxy groups (thick cir- 
cles) . The grey areas indicate the Zone of Avoidance where 
|fc| < 10°. The solid horizontal lines show the paths of the 
Arecibo beam. 



fill the beam. HVCs with column densities in ex- 
cess of the limiting value of 10^^ cm^^ could be 
detected to larger distances. The limiting col- 
umn density for clouds at large distances where 



cloud 



< c'boam is A^lim — 10 X (c?boam/rfcloud) 



For each group and galaxy, a fraction of the vol- 
ume of the surrounding sphere is scanned by the 
Arecibo beam. The number of clouds within that 
volume is calculated, taking into account the col- 
umn densities and sizes of the individual clouds. 

Table 1 lists the number of clouds that would 
have been detected in the Arecibo H i strip survey 
if a population of extragalactic HVCs existed with 
the BSTHB properties. We calculate the numbers 
of clouds differently for groups and for galaxies. 
Since BSTHB do not specify the exact radial dis- 
tribution of clouds, we tested three different ra- 
dial distribution functions to fill the volumes with 
clouds: 1) a spherical volume of radius R, 2) a 
thin spherical shell of radius R, and 3) a thick 
spherical shell with clouds distributed according 
to a Gaussian about the radius R with dispersion 
a = R/3. The latter distribution most closely 
resembles the derived distribution of Vg given in 
BSTHB. The numbers in the table are based on 



i? = 1 Mpc, the value preferred by both BSTHB 
and BB. The group halos are filled with 450 clouds, 
the number of HVCs identified by WW, excluding 
complexes A, C, M, the Outer Arm and the Mag- 
ellanic Stream. To calculate the expected number 
of clouds around galaxies, the number of clouds 
associated with each galaxy is scaled in direct pro- 
portion to the ratio of the galaxy luminosity com- 
pared to the integral luminosity of the LG. This 
leads to a median number of clouds per galaxy of 
40. 

Table 1 shows that the expected number of 
detections is essentially independent of how the 
clouds are distributed around the groups and 
galaxies. For these samples we should detect ap- 
proximately 250 clouds around galaxies and 70 
around groups if the HVC population of BSHTB is 
a common feature of nearby galaxies. Restricting 
our analysis to the compact clouds of BB, reduces 
these numbers to 39 and 9. For a uniformly filled 
spherical distribution of BSTHB clouds, the dis- 
tribution of H I masses of the; expected detections 
is shown in the lowest panel in Fig. 1. This figure 
illustrates that our analysis is sensitive to typical 
H I clouds (compare second panel of Fig. 1), and 
not only to the most massive ones. 

The robustness of the result is demonstrated in 
Table 1, where the numbers in parentheses indi- 
cate the expected number of detections if the de- 
tection threshold is increased from 5(t to 7cr. The 
average decrease is 25%. A 50% decrease occurs 
if the detection threshold is set at lOcr. This ex- 
tremely conservative threshold would still predict 
more than 100 detections. 

5. Conclusions 

The hypothesis that HVCs are primordial gas 
clouds with typical H i masses of a few x 10^ Mq 
at distances of ^ 1 Mpc from the Galaxy is not 
in agreement with observations of nearby galax- 
ies and groups. Blind H I surveys of the extra- 
galactic sky wotild have detected these clouds if 
they exist around all galaxies or galaxy groups in 
numbers equal to those suggested for the Local 
Group. These results are highly significant: the 
Arecibo H i strip survey would have detected ap- 
proximately 250 clouds around individual galaxies 
and 70 in galaxy groups. 
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Table 1 

Expected extragalactic HVCs detections 



BSTHB 



uniform shell Gaussian BB selected 



groups 70(52) 72(54) 70(52) 9(8) 14 
galaxies 248(167) 256(173) 260(177) 39(28) 347 



Note. — The number of expected HVC detections arc calcu- 
lated assuming a 5a detection threshold. The numbers in paren- 
theses indicate the expected number if the threshold is raised to 
7(7. The last column shows the number of groups and galaxies 
that are included in the analysis. For the Gaussian distributions 
a selection radius of 2 Mpc has been used for inclusion in the 
sample. 
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